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Abstract The underlying physics that generates the excitations in the global 
low-frequency, < 5.3 mHz, solar acoustic power spectrum is a well known process 
that is attributed to solar convection; However, a definitive explanation as to 
what causes excitations in the high-frequency regime, > 5.3 mHz, has yet to 



be found. Karoff and Kjeldsen| ( Astrophys. J. 678, 73-76, 2008|) concluded that 



there is a correlation between solar fiares and the global high-frequency solar 
acoustic waves. We have used the Global Oscillations Network Group (GONG) 
helioseismic data in an attempt to verify Karoff and Kjeldsen ( 2QQ8[ ) results as 



well as compare the post-fiare acoustic power spectrum to the pre-fiare acoustic 
power spectrum for 31 solar flares. Among the 31 flares analyzed, we observe 
that a decrease in acoustic power after the solar flare is just as likely as an 
increase. Furthermore, while we do observe variations in acoustic power that are 
most likely associated with the usual p-modes associated with solar convection, 
these variations do not show any signiflcant temporal association with flares. We 
flnd no evidence that consistently supports flare driven high-frequency waves. 
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1. Introduction 

Before the advent of helioseismology in the 1970s, the study of the internal struc- 
ture and dynamics of the Sun was only possible with one-dimensional theoretical 
models. While the subsequent study of solar oscillations has helped answer many 
questions pertaining to the solar interior, it has also presented new questions. 
One of these questions is whether the oscillations can be influenced by flares. 



Wolff (1972) was probably the flrst to suggest that large solar flares could excite 
free modes of oscillation of the entire Sun. Since then, much research has been 
conducted to determine if solar flares can affect the global oscillations, and in 
particular, the high-frequency waves. 

One of the flrst reports of the effects of solar flares on solar oscillations was 



published by Haber, Toomre, and Hill (1988). Dopplergrams of intermediate- 



degree p-modes were used to observe the change in acoustic power before and 
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after a major solar flare event. They concluded that the flare caused an average 
increase in power the day after the flare. Shortly after, [Braun and Duvall (1990) 
analyzed the interaction between the p-mode oscillations and an X-class flare 
in NOAA (National Oceanic and Atmospheric Administration) active region 
5359 using radial velocity observations obtained from the Kitt Peak vacuum 
telescope. Their results indicated no evidence of an increase in the oscillation 
amplitude correlated with the flare event. Foglizzo ( [l998 ) studied the corre- 
lations between low-degree modes to determine if they could be excited by 
mechanisms other than granulation. It was concluded that X-class solar flares 
most likely do not excite low-degree oscillations due to inconsistencies in energy 
conservation. Ambastha and Antia (2006) search for a correlation between the 



power of low-degree modes observed in velocity and flare energy release was 
generally inconclusive. Gavryusev and Gavryuseva^ (1999) analysis found a high 
correlation between temporal changes in low-degree spectral power and coronal 
mass ejections, which suggests that there exists a relationship between p-mode 
power and solar activity. Although, it is unclear whether thier finding can be 
extended to solar flares. Donea and Lindsey (2005) found that a helioseismic 
holography analysis of the seismic emission from two X-class solar flares indicates 
that the hard X-ray emission produced during solar flares are associated with 
acoustic signatures. Kumar and Ravindra (2006) reported an enhancement in 
power maps at frequencies from 5 mHz to 6.5 mHz of a solar active region in 
which a large X-class flare was produced. Finally, Kosovichev has observationally 



shown that flares can excite local waves at the location of the flare ( [Kosovichev 
and Zharkova||1998| |Kosovichev||2006[ |Kosovichev||201lD . 



Thus, while it is established that large flares can generate local acoustic waves 
in the solar atmosphere, there is still conflicting evidence for the role of flares in 
affecting global oscillations. In particular, "Karoff and Kjeldsen (2008^), hereafter 
KK, used solar disk-integrated intensity data from the Variability of Solar IRra- 
diance and Gravity (VIRGO) instrument, on-board the Solar and Heliospheric 
Observatory (SOHO) satellite, and solar X-ray flux data from the Geostationary 
Operational Environmental Satellites (GOES) to study the effects of flares on the 
global oscillations. They found a high correlation between the power in the high- 
frequency acoustic power spectrum and variations in the solar X-ray flux. Kumar 
and collaborators have found evidence supporting the theory of flare driven high- 
frequency acoustic oscillations ( Kumar and Venkatakrishnan|2009[ [Kumar et al 



2010). Using a wavelet analysis of Michelson Doppler Imager (MDI) and Global 



Oscillation at Low Frequency (GOLF) data, they found enhancements of high- 
frequency low-degree power associated with three powerful flares in the MDI 
data, but only marginally in the GOLF data. Both of these instruments are 
also on SOHO. In addition, Kumar et al (2011) reported a significant post- 
flare enhancement in the high-frequency global velocity oscillations for both 
GOLF and the Global Oscillations Network Group (GONG) velocity data of a 
X-class solar flare in NOAA active region 10930. Most recently, [Chakraborty[ 



et al (2011) were able to reproduce the KK results using acoustic data obtained 



by the Sun PhotoMeter (SPM) instrument on SOHO and GOES X-ray data 
obtained by the Solar X-ray Imager. However, the results they obtained from 
correlating the GOES X-ray data with the acoustic data produced by GONG 
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instruments showed a weaker correlation between solar flares and acoustic power 
at frequencies higher than the solar acoustic cut-off frequency. 

In this paper we present an analysis using GONG and GOES data to deter- 
mine if global helioseismic signatures are affected by major solar flare activity. 
The analysis was conducted using two independent methods. The first method 
was a comparative analysis between pre-flare and post-flare acoustic power spec- 
tra for a few selected major solar flares. The objective of this portion of the 
analysis was to see if there were increases in the power spectra directly after the 
flare event as Wolff's (1972) hypothesis suggested. The second method was an 
attempt to reproduce the KK results by constructing frequency vs time plots for 
four different ranges of spherical harmonic degree and comparing those results 
to temporal variations in the solar X-ray flux. Given that space-based detectors 
are constantly being bombarded by energetic solar particles from events like 
solar flares and coronal mass ejections, it is possible that the KK data set may 
have been affected by such particle events, which could influence the data. As 
GONG data are ground-based solar observations and can not be affected by 
such interference, the GONG observations are an ideal data set to use to test 
the results obtained by Karoff and Kjeldsen (2008). 



2. Data and Analysis Methods 

In this study we examine approximately 14 years, from Dec. 9, 1995 to April 
29, 2010, of GONG global helioseismic observations. The GONG instruments 
generate polarized disk intensity images of the Sun, taken every 60 seconds, 
as their primary data product. After reductions and calibrations, the intensity 
images are converted to Doppler images. Each Doppler image from each of the 
six GONG sites is decomposed into a set of spherical harmonic amplitudes or 



coefficients ranging in degree I from to 200 (Hill et al., 1996). The coefficients 
from simultaneous measurements are merged into a set of coherent time series 
for each individual value of i and |m|, where m is the azimuthal degree, with 
< |m| < I. The spherical harmonic time series data show the temporal behavior 
of global solar p- modes. 

The time series are converted into power spectra by taking the FFT of the 
corresponding spherical harmonic coefficient time series. This produces spectra 
as a function of m and v for each ^, with —i<m<i. This step also separates 
the modes in frequency u with identical {£^m) but different n, where n is the 
radial order of the mode. The azimuthal order m yields information about the 
phase and spatial pattern of a p-mode on the solar surface. However due to solar 
differential rotation, multiplet p-modes with identical (n, i) but different m that 
would have identical frequencies if the solar surface rotated as a solid body, 
acquire non-degenerate frequencies with respect to azimuthal order m. Note 
that, since we cannot observe the entire spherical surface of the Sun, a power 
spectrum computed for a specific £ also contains power from several adjacent 
values of ^, producing additional features in the spectrum. These artifacts are 



known as spatial leaks (Howe and Thompson, 1998). 
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In order to decrease the noise in the spectra generated by the stochastic 
excitation of the modes, we next shift the spectrum for each n) multiplet by 
a frequency 5u as a function of m by the amount Su = mP(£ ,^), where P(£^n) is a 
Legendre polynomial. This removes the frequency shift arising from differential 
rotation and other effects. Next, the spectrum is averaged over m, and then over 
four ranges of as follows: £i: < £ < 50; £2: 51 < £ < 100; £3: 101 < ^ < 150; 
and £4: 151 < £ < 200. These m- and averaged spectra are computed over a 
variety of time spans, as described below. 

We used the soft channel (1-8 A) flux measured by the X-ray instruments 
on the GOES satellites to compare the occurrence and strength of major solar 
flares to increases in solar acoustic power. This time series is averaged over the 
same time intervals as the oscillation data for comparisons. Rapid changes in 
the X-ray flux are used to identify the times of flare maximum. 



2.1. Pre-flare/Post-flare Comparison Analysis 



We compared oscillation power spectra before and after individual large flares to 
test Wolff's hypothesis ( 1972 ) that flares can stimulate the normal modes of the 
Sun. We selected 31 flares that were of M5-class or higher and within 45° of the 
center of the solar disk. For each flare, we calculated pre-event and post-event 
average power spectra of one day duration, divided at the time of the maximum 
intensity of the flare. The power ratio of the post-flare to pre-flare power spectra 
was then calculated as a function of u in the four ranges of £ described above. 
As a control, the same procedure was performed for points in time when major 
solar flares, M5-class or higher, were not present. 



2.2. Production of power as a function of frequency and time 



For confirming the KK results, we computed power as a function of frequency 
and time in a manner similar to the method described in Karoff and Kjeldsen] 
( |2008) , but with a few differences. Consecutive 7.5-day m-averaged, ^-averaged 
FFT power spectra were calculated as explained above. We chose 7.5-day time 
strings so that the frequency resolution of our analysis was in agreement with the 
frequency resolution of the KK analysis. Each time string was displaced 18 hours 
from the start time of the preceding time string and the power spectra were then 
vertically stacked to produce the plots. KK's analysis employed smoothing of the 
power spectra. No smoothing was employed in our analysis so as to preserve the 
full rich information content of the data especially for the high frequencies of 
interest here. This procedure was performed for the four different ranges of £. 
In addition, we further averaged these data over frequency and over time to 
examine the overall variations of power in time and frequency, respectively. 

For consistency, the X-ray data were processed in the same manner as the 
power spectra. The X-ray time series was averaged in 7.5 day time strings in 
which the starting time of each interval was shifted by 18 hours from the starting 
time of the previous time string. No smoothing was applied to the X-ray data. 
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3. Results 



3.1. Pre-flare/Post-flare Comparison 

Given the number of flares used for this portion of our analysis, we present plots 
for only three X class solar flares selected for our comparison analysis. Results 
obtained for the entire set of flares are summarized in Table 1 and graphically 
displayed in Figures |5]and[6l In Figure [ij for the X5.7 flare of July 14, 2000 (the 
Bastille Day flare), we see a general increase in acoustic power after the flare 
for all four ranges of ^ used for this analysis, and for frequencies above about 2 
mHz. For this case, Wolff's hypothesis could be considered to be supported. 

However, in Figures |2] and |3] the plots generally show a decrease in power, for 
all ranges of I used, after the occurrence of the flares. Both of these flares, an X6.2 
on December 13, 2001, and the Halloween X17.2 event of October 28, 2003, were 
substantially more powerful than the Bastille Day flare of Figure 1, but show a 
relative decrease in power rather than an increase. Of the 31 solar flare events 
selected for our analysis, as shown in Table 1, 16 (52%) have an average power 
ratio above 1 in the ii range; 13 (42%) have an average power ratio above 1 
in the ^2 and ^3 ranges, and 14 (45%) have an average power ratio above 1 in 
the ^4 range. The plots in Figure [4] are for a control quiet time, and show an 
overall decrease in power similar to that in Figures [2] and [3] From these results, 
we can conclude that flare-related variations are probably no different from the 
variations during quiet times that arise from the stochastic driving of the modes. 
Furthermore, we calculated the correlation for the data presented in Figures [5] 
and [6] using both Spearman and Pearson correlation statistics for all four of the 
^ ranges plotted. The Spearman rank correlation yielded correlation values less 
than 60% confidence and the Pearson correlation values strongly indicate no 
correlation. 

There is a clear increase in the variability of the power ratio above the acoustic 
cutoff frequency around 5.3 mHz for the lowest range ^1 (0 < ^ < 50). This 
excess variability is observed for almost all of the 1-day spectral ratios, both for 
flares and for the quiet Sun. In addition, the high-frequency power ratio in the 
^1 range is always much more variable than at the higher degrees, as shown in 
Figure |5| While this suggests that there may be a difference in the excitation and 
damping of the lowest degree modes, there are also far fewer values of m available 
for averaging at the lower degrees. Thus, the higher variation with respect to 
frequency for the lowest range ^1 may simply be the result of incomplete noise 
cancellation. However, the variability in the power ratios for all three of the 
higher ranges of ^ do not show an increasing reduction that would be expected 
if the effect is mainly due to m averaging. 
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3.2. Power as a function of frequency and time 



Figure [7| shows the average p-mode power as function of time and frequency in 
our four ranges of spherical harmonic degree i. The prominent red/yehow bands 
are indicative of the usual p modes, which have typical frequencies between 2.2 
and 5 mHz. The left side of Figure [7| shows the GOES soft X-ray flux as a 
function of time, averaged over the same 7.5 day periods offset by 18 hours that 
were used for computing the power spectra. Note that the four color plots have 
different scales. 

A comparison of the location of the peaks in the soft X-ray flux with the 
localized variations in the oscillation power spectra does not show any obvious 
correlation between the X-ray flux and enhancements of oscillation power. The 
lack of correlation can also be seen in Figure [8) where the data in Figure [7| 
have been averaged over all frequencies, producing the temporal evolution of the 
oscillation power in our four ranges of i. 

Both Figures [7| and |8] show a number of interesting features. The large drop 
in power in the ii range in early 2000 is very likely a data processing artifact, 
as is the similar drop evident in Figure [7] in the ii and £2 ranges in late 2004. 
A significant number of outliers were discarded from the power spectra at these 
times. 

There is a sharp local increase in power in all ranges of i in April 2001, which 
follows a series of three major flares within one month. These flares were an XI. 7 
on March 29, 2001; a X2.3 on April 10, 2001, and an M7.8 on April 26, 2001. 
It is possible that this sequence may have indeed amplified the strength of the 
oscillations. Finally, Figure [8] shows that there is an abrupt and then persistent 
change in the overall power level of the oscillations in the is and £4 ranges in 
mid 2005. The change is also apparent in Figure [7| as a change in the width in 
frequency of the p-mode band. 

The origin of this change is most likely due to a revision of the method used to 
reconcile the angular orientation of the solar images around the GONG network. 
Since it is impossible to perfectly align the instruments north-south on the earth, 
the apparent position of the solar rotation axis varies by a few tenths of a degree 
between the six GONG sites, and this discrepancy changes over the course of the 
year. An angular change in the rotation axis orientation redistributes the power 
attributed to a given oscillation mode with specific values of i and m into other 
values of £ and m, and this introduces errors into measurements of the frequency 
and amplitude of the modes ('Kennedy 1997 Kennedy and Williams|1998 ). These 
errors are larger at higher values of ^, and also increase with an increase in the 
rotation axis discrepancy, so considerable effort has been made to reduce the 
site-to-site angular discrepancies in the GONG data and to determine the true 
orientation of the Suns rotation axis. It is estimated that these efforts have 
decreased the systematic error in the axis orientation by an order of magnitude, 
down to 0.02°. In 2005, the GONG Project introduced a scheduled drift scan 
observation at local noon around the network that substantially improved the 
accuracy of the orientation measurement. This is the source of the amplitude 
drop in 2005 seen in Figure [Sj which shows no drop for the lowest range of ^, 
and an increase in the drop as i increases. This is consistent with the known 
properties of the effect. 
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Figure [9] shows the corresponding time averaged power spectrum as a function 
of frequency for the four ranges in £. We see that the acoustic power at frequen- 
cies higher than the acoustic cut-off frequency around 5.3 mHz are greatest in 
the £i range (-^ = to = 50) (top left panel). This suggests that the solar 
processes that excite the high-frequency oscillations do so mainly at low and 
intermediate degrees. The time-averaged spectra for the is and £4 ranges in 
Figure [9] display some notch-like jumps, which are also the result of the change 
in the determination of the solar axis orientation as discussed above. As before, 
Figure [9] shows that the magnitude of the change increases with the range of £ 
and, in addition, that the power redistribution occurs at frequencies that also 
depend on the range of £. 



4. Discussion 



In this study, we have searched for evidence that large and energetic solar flares 
can change the amplitude of the global oscillations of the Sun, as postulated by 



Wolffl(|l972) 
)8l). 



and apparently observed in SOHO data by Karoff and Kjeldsen 

(|2008|). We find no evidence for this process in the GONG data. Ratios of m- 
and averaged spectra obtained before and after 31 flares of class M5 or higher 
essentially show that a decrease in the mode amplitude is just as likely as an 
increase, and images of the mode power as a function of time and frequency 
do not show any significant variations in time that can be associated with 
flares. Our results do show variations that are most likely the result of the 
stochastic excitation of the solar p modes. There are some puzzling facets of 
this variation, such as the relatively large amplitude variations at low spherical 
harmonic degrees below ^ = 50, compared to the variations at higher degrees. 
The variance of the mode amplitude is also observed to be higher at all £ above 
the acoustic cutoff frequency of 5.3 mHz. This may be due to the lack of trapping 
of waves with frequencies higher than the acoustic cutoff, which could allow the 
rapidly evolving solar atmosphere to have a larger influence on the coherence of 
the waves compared to the relatively quiescent interior. 



With the exception of [Kumar et al. (2011), other studies have not detected 
flare-related changes in the global mode amplitudes as observed by GONG 
(Chakraborty et al. 2011 Ambastha and Antia 2006), but studies using data 
from instruments on SOHO generally have detected such changes (Karoff and 



Kjeldsen! [20081 [Kumar and Venkatakrish nan| [ 20091 [Kumar et aJ] [2010|[Chakrabo| rty| 



et 



2011). There are at least five possible explanations for this discrepancy. 



One is that the noise level of the GONG observations is substantially higher 
than the noise in the SOHO instruments, but there have been many comparisons 



between the data sets since 1995 (see, e.g., Komm et al, 1998), and there is no 
indication that this is the case. A second possibility is that all of the instru- 
ments used by various studies of this question have different distributions of the 
spatial leaks discussed in Sec. 2. These leaks may somehow mask the amplitude 
variations associated with flares in the GONG data. The third possibility is 
that the flare effects are restricted to only the very low-degree modes, and that 
averaging over a range in £ of to 50 has masked the effects. On the other 
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hand, it would be surprising if a very localized impulse such as flare would 
only affect very low spatial frequencies, since the spectrum of an impulse in 
the Fourier domain is essentially white, with all spatial frequencies present. The 
fourth possibility is that the SOHO instruments are directly affected by the 
flares since they are located in space and exposed to the excess particle emission 
generated by these events. Ground-based instruments are shielded from this by 
the Earth's atmosphere. It is possible that the apparent changes in the oscillation 
amplitudes are due to increased noise when charged particles hit the space-based 
detectors. A test of this would be to search for flare effects using BiSON data, 
which is another ground-based network. The fifth possibility may be related 
to the way in which the data have been reduced and analyzed. An analysis of 
MDI, GOLF, and VIRGO data using the methodology of this analysis would 
be an ideal way to test for discrepancies due to data processing. Although the 
aforementioned analyses are outside the scope of this effort, it is possible that 
they will be done for a future paper. 

The results also raise some additional questions about the interaction of the 
flares and the global oscillations. The flare associated with the largest increase 
in amplitude in Table 1 is the M5.9 flare on October 5, 2002. The power of the 
oscillations basically doubled for this event at all degrees, even though it was 
a relatively weak event. Why? Was this simply a coincidence between the flare 
and a strong random fluctuation at all values of or was there some special 
aspect of the event, such as its location on the solar surface? Does the location 
or time of an event play a role in how effective a flare may or may not be in 
modifying the oscillations? Finally, does a sequence of flares, such as occurred 
in March and April of 2001, have a cumulative effect on the modes? If so, how 
often and how strong do the flares have to be in order to cause such an effect? 
With the advent of asteroseismology, we may be able to address these questions 
using statistical studies on data sets from Kepler and the Stellar Observations 
Network Group (SONG). 
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Figure 1. The power ratio of m-averaged, -^-averaged temporal power spectra before and after 
the X5.7 flare that occurred on July 14, 2000. The ratios are in the sense of post-flare power 
divided by pre-flare power, so a ratio greater than 1 indicates a power enhancement. 
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Figure 2. As for Figure 1, but for the X6.2 flare of Dec. 13, 2001 
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Figure 3. As for Figure 1, but for the X17.2 flare of Oct. 28, 2003 
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Figure 4. As for Figure 1, but for an arbitrary quiet period on April 3, 2004. 
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Figure 5. The average value of the power ratio above the acoustic cutoff frequency (> 5.3 
mHz), as seen in columns 4 - 7 in Table 1, plotted as a function of solar X-ray intensity. The 
error bars are the corresponding standard deviations in columns 8 - 11 in Table 1. 
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Figure 6. Percentage of values greater than 
mHz), as seen in columns 12 - 15 in Table 1, 
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Figure 7. The color images show the averaged oscillation power spectrum as a function of 
year and frequency v for the four ranges of i. The color scale shows the magnitude of the 
power and is logarithmically scaled. Note that the color scales are different for each of the 
images. The red/yellow vertical bars starting around 2 mHz and ending around 4.5 mHz in 
the color images are the p-mode oscillations. The sharp decrease in power seen in the top left 
panel around early 2000 and late 2004 as well as in the top right panel around late 2004 is 
due to outlier removal. The two narrow panels on the left of the figure show the average soft 
X-ray flux as a function of time. The large horizontal spikes correspond to major solar flare 
occurrences. 
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Figure 8. The temporal variation of p-mode power averaged over all frequencies. The X-ray 
flux is shown the in bottom left and right panels. There does not appear to be any consistent 
relationship between X-ray flux and acoustic power. 
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